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While vector modes are usually ignored in cosmology since they are not produced during inflation 
they are inevitably produced from the interaction of density fluctuations of differing wavelengths. 
This effect may be calculated via a second-order perturbative expansion. We investigate this effect 
during the radiation era. We discuss the generation mechanism by investigating two scalar modes 
interacting, and we calculate the power of vector modes generated by a power-law spectrum of 
density perturbations on all scales. 



I. INTRODUCTION 



• High-precision data from observations of the Cosmic Microwave Background (CMB) and Large Scale Structure 
' (LSS) [2] provide strong evidence for a nearly spatially flat universe with a primordial spectrum of adiabatic, Gaussian 
I and nearly scale-invariant density perturbations. The standard cosmological model provides a remarkable theoretical 



basis for these and other observed features of our universe. Perturbations generated from inflation give a nearly 
^ scale-invariant spectrum of scalar (density) perturbations, and tensor (gravitational wave) modes, with amplitudes 
I which are typically within a few orders-of-magnitude of each other. However, within this paradigm the amplitude of 
, any corresponding vector modes is zero since a scalar field can't support vector modes at linear order ; even if they 
were generated during inflation, vectors decay rapidly after they leave the Hubble radius during inflation, whereas 
I— 1^ scalars and tensors are typically frozen on super-Hubble scales, and only decay when they re-enter the horizon. A 
i-S^ . generic prediction of inflation is therefore no vector modes. 

There is an important caveat to this argument. Vector modes are generated via the non- linear interaction of scalar 
Q . (and tensor) perturbations of differing wavelength, and therefore inflation must generically predict a spectrum of 
' vector modes, but at second-order in a perturbative expansion. Indeed, observations of the scalar spectrum requires 
c/2 , this to be so independently of whether inflation is the correct model of the early universe or not. We shall consider 
, ^ , ■ the generation of vectors from scalars in some detail. The analo gou s process of gravitational wave generation by 
■ scalar-scalar interaction has been investigated 0, S S, 11 B [Mllil, [l3, [lE El, 111, d, [13, E Ei , and the work 
PsJ • presented here is closely related to these studies. 

K*" ' We shall principally expand on work of MoUerach et al. [20] , who considered the effect of secondary vectors on the 
CMB. We shall investigate the generation method of vectors from two scalar modes (we show that unlike gravitational 
waves, vectors can't be generated by a single mode), and we shall calculate the power spectra of the vector part of 
the metric in the radiation era. Our aim here is to principally discuss how vectors are generated, and get an overall 
estimate as to the magnitude and distribution of vectors at the end of the radiation era. How the results here relate 
to observables in the CMB, and the spectrum of vectors today, is left for future work. 

Other people have discussed second-order vector modes before. Most recently, Mena et al. [2l| considered second- 
[ order vector modes in a collapsing universe. Matarrese et al. p2| discussed the generation of primordial magnetic 
fields from density perturbations (although see [23j). Various other work has discussed vector mode generation on a 
more formal level [ll [H, [13, [H, [ll, [13, [21] . 

There are a variety of other mechanisms which predict vector modes, all of which must happen after inflation, and 



o 



X • 

5_( ■ usually predict a spectrum of modes on small scales. Such sources include cosmic strings [29|, topological defects 3C|, 
5^ ' fine-tuned anisotropics in coUisionless neutrinos [sl*! and the presence of an primordial magnetic field [13, [H, [H, 35 1 . 

The generated vector modes are highly non-Gaussian. Second-order inflationary vector modes considered here are 

also non-Gaussian and have a x^-distribution. 

Vector modes are likely to play a more prominent role in cosmology in the coming years, through their contributions 

to the CMB [13, [25, 26, 36, 37, 38, 39, 40], which will impact on the B-mode polarization and could be the dominant 

contribution when compared with the second order gravitational waves |13] (assuming one is able to subtract the 
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lensing signal). They are also crucial in magneto-genesis [H, [H, [H H [H, SH, US, H H SI, IH, 113, El, Ell as 
the magnetic field vector has a dominant vector part. It is therefore appropriate to now consider in more detail the 
spectrum of vector modes which we know must exist by virtue of second-order effects. 

The paper is organised as follows. In Sec. II we consider the formalism for investigating the generation of vectors 
from scalars. Then, In Sec. Ill we discuss the power spectra of vectors in the radiation era for power-law scalar 
modes. We also investigate the generation mechanism of vectors by considering the interaction of two distinct scalar 
modes. Finally, we conclude in Sec. IV. 



II. GENERATION OF VECTOR MODES FROM DENSITY PERTURBATIONS 



We shall consider perturbations of a flat Robertson- Walker background up to second order. The metric is decom- 
posed as 

= 5m"^ + "^^^^Sm"^ + "^^^^Sm"^' (1) 

where Greek indices run from ... 3 and Latin indices run from 1 ... 3. We are only investigating the second order 
vector modes sourced by the first order scalar perturbation; therefore, we have scalar degrees of freedom at first order, 
6^^^ g^y, while the second order perturbations, (S'-^-'g^jy, are pure vector modes, and second-order equations are therefore 
meaningful when projected out accordingly. Our metric in terms of conformal time with a longitudinal gauge chosen, 
is given as 

5oo = -a'(l + 2$(i)), 1^2^(2)^ = (l - 2$(i)) (2) 

where <^'^^^ is the first order Bardeen potential, and sf'^ describes the gauge-invariant [2l| second order vector modes, 

(2") 

SO that d'^Sl =0. As there is no ambiguity in what follows we shall drop the order superscripts and just write <I> 
and Si. 



A. Density perturbations at linear order 

We shall consider the generation of vectors from scalars during the radiation era below, but for now we shall 
consider the situation where we have a linear constant equation of state, = w where w = p/ p and is the sound 
speed. We shall assume that the first order matter perturbations are adiabatic, i.e., that pressure perturbations obey 
5^^)p = c^(5^^^p. Then the first order equation of motion for the Bardeen potential in Fourier space is [50l | 

<i>"(fc, 77) + m{l + c2)$'(fc, r;) + clk^Hk. V) = 0, (3) 

where a prime denotes differentiation with respect to conformal time rj. 

In the radiation era, the scale factor, the conformal Hubble rate and energy density evolve as a(r]) = oo^i ^ = 

— = - and p (X 1]^^, and the general solution to ([3]), with = -i is 
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We shall ignore the decaying mode - that is terms with a Br{k) coefficient. 

Assuming that the fiuctuations are Gaussian, we may introduce Gaussian random variables, E, with unit variance 
and the property 



(5) 



We can then separate the length and directional dependence of functions of k and write $(^,77) — ^{k,ri)E{k) and 
Arik) = Arik)E{k). 

The power spectrum for the first order scalar perturbation can be defined through 



($*(fci,r?)$(fc2,r/)) = ^S^ki ~ k2)V^{k,7j). 



(6) 



At early times during the radiation era the power spectrum becomes 

4867r 

Relating the Bardeen potential to the comoving curvature perturbation at early times gives us 



r^ik) ~ Arikfj^. (7) 



« ^A^(fc), (8) 

where is primordial power spectrum for the curvature perturbation TZ. Current observations show w 2.4 x 10~^ 
at a scale kcMB — 0.002Mpc"^, and is almost independent of wavenumber on these scales ^ . 



Second-order vector modes 



The vector perturbations at linear order satisfy an evolution equation and a momentum constraint equation which 
can be found by calculating the i — j and — i parts of the Einstein Field Equations (EFE's) respectively [2l|. In 
the case of a perfect fluid (at first order only) there is no source in the vector evolution equation and it admits 
solutions proportional to 1/a^ . The momentum constraint equation relates the vector perturbation to the 3- velocity 
perturbation, which in the perfect fluid case would be zero. However, the respective equations at second order differ 
significantly. Firstly, as we will see, the evolution equation is sourced, allowing for the generation of vector modes. 
Secondly, the momentum constraint no longer excludes the existence of vector perturbations, provided we are only 
considering a perfect fluid up to first order (see [2ll|). 

We calculate the evolution equation for Si in the usual manner, by expanding the EFE's up to second order, keeping 
terms quadratic in the first order quantities. We start with the trace reversed EFE's 
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The second order space-space part of the Ricci tensor can be written as 



where we have 



which contains only the second order term; and 
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which are the quadratic first order scalar perturbations. 

The second order trace reversed space-space part of the energy momentum tensor is 

<5(2)T,, = a" (1 - w) $<5(iV7., + (1 + w)p {d,V(^)) (a,i;(i)) 
It is obvious that we require the following zeroth and first order equations. 



(13) 



= -nGa^p, 
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Note that the terms with as a coefficient will not play a role since the terms are eliminated by the operator 
Vj-™ defined below. Due to the limited quantities we keep in our metric the tensorial equations we calculate are only 
valid for vector modes, and so these must be projected out. 



We define the Fourier transform of the vector perturbation as 

S^{x, ,7) = ^^-L^ j d?k [S{k, ri)e,{k) + 5(fc, 'n)e,{k)\ e*-, (14) 

where the two orthonormal basis vectors e and e are orthogonal to k. We shall use the operator V,'™ to extract out 
the divergenceless vector from a rank-2 tensor 

V,'" = --^ ( <fk' I ^3^' yi [eYA;')e™(fc') + e,(fc')e"(fc')] e*'=' (^-=='). (15) 
(27r)'^ 7 7 

This operator will remove any rank-2 tensor which is constructed from derivatives of a scalar potential. Specifically, 
any second order scalar perturbations will be removed and therefore such perturbations have been neglected in our 
analysis. The action of this operator is to produce a rank-1 vector which is a pure vector mode. Further details 
concerning the operator can be found in Appendix [XI 

We can then obtain, from the i ~ j component of the EFE's, the evolution equation for the second order vector 
perturbations 

K"" {d(iS[n) + ^nd^iSm)) = 2V^5],™, (16) 

where the source term is given by 

E,^ = -4$ (9,9™$) - ^^^-^ (9,$) (a„,$) + 1 [(9,$') (a„,$') + 2H (9,$) (a„,$')] ■ (17) 

For either polarisation, the evolution equation in Fourier space for the vector mode becomes 

S'{k,'q) + 2HS{k,Tj) = Y.{k,Tj), (18) 
where the source term S(fc, 77) is an appropriate convolution over the quadratic first-order quantities, 

10 + 6w 
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3(1 + w)n^{r)) ^ ' '" ^ ■" ' ' 3(1 + w)H{r)) 
The general solution for S (fe, 77) can be written as 
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We have set the initial conditions for the vector mode to zero at 7; = 770. 

1. Power spectrum 

The power spectrum of the induced vector mode is defined as 

27r2 



(19) 



5(fc,77) = ^--/ d77a^(r7)S(fc,77). (20) 



{S* (fci, 7;) S (fc2, 77)) = (fei - fe2) Vvik, 77). (21) 

Substituting the solution ([^0]) into (PT|) we find 

{S* {k,,ij)S{k2,v)) - f df,2 r dm a'(7?i)a2(772)(S*(fei,^i)S(fc2,7?2)) 

16 [A;i"e"(fci)] [fc2V(fc2)] f, , 2,^,2,-, 

' ari2 / 0771 a [rjija (7^2) 



(27r)3a4(77) ^ 'V,o " ' '''' (22) 

X j d^k[j d'k'^ (fcUfciJ(fc^,fc^,)S(fci,|fci -fell, 771) S(fc^,|fc2-fci|, 772) 

X (E*{k[)E*{ki - k[)E{k'2)E{k2 - k'^] 



where 

However, Wick's theorem tells us that 

[E*{k\)E*{ki - k[)E{k'2)E{k2 - k'^)) - {E*{k[)E*{ki - k[)) l^E{k'^)E{k2 - k'^) 

'e* (fci)^(fc^)) (e* (fci - fc'i)^(fc2 - fei) ) (23) 
'E*{k[)E{k2 - k'^)) (E*{ki - k[)E{k'2) 



Therefore, the power spectrum of the induced vector mode is 
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Vv{k,v) = ^^5^4(^) J J dfj, a'im)a\f„) J d^k'ik^K) [e\k)K] [e^ {k)^] ^^^^ 

X S (fc', \k - fc'l, 771) [(fc'fcO S (fc', |fc - fc'l, 772) + {k'ki - S (|fc - k'\,k', fi2)] . 
In order to compute the integrals over Fourier space, we first introduce the diniensionless variables u and w, where 

k' 

V — — and u = \'l + v'^~2v cos 6. 
k 

If we rewrite equation (|24p using spherical coordinates in Fourier space, we can carry out the azimuthal integral 
trivially. Using the two new variables, the power spectrum then becomes 

^8 rv rn roo /-u+l 

Vv{k,7l)^ 1. 4 4/ W dm dma\m)a^f,2) dv / du{uv){v^ + l-u')[{u^^l-v^f~Av^] 

iOTT a [rj) j^^ Jq J\v-i\ (25) 

X S {kv, ku, fji) { (u^ — 1 — w^)S {kv, ku, 7)2) + ('"'^ + 1 — w^)S (fcu, /cw, 7)2)} ■ 
The power spectrum can now be calculated once the power spectra (initial conditions) for the scalar modes are chosen. 

III. VECTOR MODE POWER SPECTRA 

We shall now investigate the power spectrum of the induced vector modes during the radiation era. 

After substituting for the first order solution for $ for the radiation era, the power spectrum then becomes 

(243)2 /-ao 

Vy^k,!]) — — — -J / dv du V<s>(ku)V,i>(kv)J^(u,v,x), (26) 

4(fc77)4 7o 



where 
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T(u, V, x) = (t;^ + 1 " u^) liu^ ~ 1 - v^f - Av^ 
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and X is another dimcnsionless variable defined as a; = ^77, and xq — krjQ. We have defined the functions 

U^) = E E (""^ + '^«) -^1^ (28) 

m 

with the coefficients a^, 4>n and Af^^ defined as 

'^""2(0; — " x/3 V ul^va ; ' 
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Here we have defined the matricies 1, 0, a and b as 
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As we have four integrals to carry out it is useful to calculate them analytically where possible. We can do this for 
the x-integrals to get 
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For the radiation era we assume that all modes are well outside the horizon when the interaction begins and therefore 
can set xq = 0. 



A. Interaction of scalar modes 



Before calculating the power spectrum for the case of power law scalar modes, it is useful to investigate how the 
vector modes are generated from individual scalar modes. It has been shown that a single scalar mode with an 
isotropic distribution will induce second-order gravitational waves This is not the case with vector modes: scalar 
modes of differing wavelengths need to interact to generate vector modes, as we shall see. To investigate this we 
choose then a scalar power spectrum of the form 

P$(fc) = ^A^A^T^ikcMB) {^[ln(fci/fc)] + S[\n{k2/k)]} , (34) 

where A is the mean amplitude of each wavenumber, ki, relative to the observed amplitude of the primordial power 
spectrum, A'^{kcAiB), at wavenumber kcMB 3> ki. We assume for simplicity that they both have the same amplitude. 
Carrying out the u and v integrals, we then find that the vector power spectra becomes, in terms of Vi = ki/k, 

rvik,v) = 262UA*A^'^ [Tivi,V2,x) + Tiv2,vux)] (35) 

X K 

provided ui + 1 > 7J2 > |wi — 1|, and is zero otherwise. Therefore modes are induced for all wavenumbers k such that 
fci + fc2 > A: > |fci — fc2|, and are scattered into angles such that mi = W2, and 1*2 = fi (a further requirement from 
carrying out the integration). In the case where only one input mode is present, so V2 = wi, this inequality becomes 
ki > fc/2 while we also have vi — ui, as in p^ : from Eq. (j27p with u = v, however, we see that vanishes in this 
case. Therefore, we can see that vector modes cannot be induced by a single scalar degree of freedom. The physical 
reason for this is because vector modes are associated with rotational degrees of freedom. A consequence of ki — k2 is 
that 9i — ±02, i.e., the input modes only have momentum along the same axis in Fourier space. Consequently, there 
is no angular momentum generated, and hence no vectors. 

Provided that fci 7^ fc2, we can have vectors induced over the appropriate range of wavelengths. Closely separated 
scalar modes will produce a much broader spectrum of vector modes while modes of vastly differing wavelengths will 
produce a very narrow range of vectors, with wavenumbers close to the largest input wavenumber. Note that as the 
generated wavenumbers are restricted from above and below, we can't expect any noise on large scales, as is the case 
for gravitational waves. This is also evidenced by the fact that one input mode can't produce any vectors - there 
would be nothing to set the long wavelength cutoff in that case. 

In Fig.[l]we show the induced vector modes as a function of x for various vi with V2 — ^ (so we require < wi < 2), 
i.e., the evolution of modes of wavenumber k — k2. While the generated mode is outside the Hubble radius, there is 
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FIG. 1: The power spectrum of vector modes induced by two interacting scalar modes. Although maximum power is generated 
in the scenario k\/k2 ~ 1 shortly after the Hubble radius is crossed, at late times scalar interactions with vastly differing 
wavelengths produce more power once the long- wavelength mode enters the Ifubble-radius. 



power-law growth, with Vv ~ rf' . When ^2 enters the Hubble radius, the principle generation of vector modes stops 
shortly thereafter, and the induced modes start to decay as ?7~^. This continues until the longer wavelength mode 
enters the Hubble radius at kir] = 1 x ~ l/wi- This then generates a further burst of vector modes, which we 
can see by progressively more pronounced knees, as vi — *■ 0, in the curves at late times. For the case when vi — 1.5, 
on the other hand, we see some confusion as the modes enter the horizon more-or-less together, before decaying as 
normal. 

Thus we see that vector modes are only efficiently generated when at least one of the scalar modes is entering its 
Hubble radius - provided another scalar mode exists to help seed the vector mode. This explains why we have knees 
in the evolution of the generated vector modes, since two interacting scalars enter the Hubble radius at different 
times. The power generated into vectors as each mode enters depends on the relative ratio fci/fc2. Modes of similar 
wavelength generate more overall power, because they are entering at the same time; modes that are widely separated 
in wavenumber don't generate as much overall power but produce more pronounced knees instead. 



Let us now investigate the spectrum of vector modes from power-law scalar modes. To do this, we assume that the 
input power spectrum is: 



B. Power law scalar modes 




(36) 



where the index Us tehs us the tilt of the spectrum relative to scale-invariance, ris = 1, and kcMB is a pivot scale for 
the power spectrum p_] . The induced vector modes are then given by 

where J^n^ (^) is defined as 

Tnjx)= / dv du {uvY^~^T{u,v,x). (38) 

Jo J\v-1\ 

We integrate this numerically, and show the results in Fig. [2] for the case Ug ~ 1. The tilt of the scalar power 
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FIG. 2: The power spectrum of vector modes induced by scale-invariant scalar modes. Scalar modes outside the Hubble radius 
interact to give power-law growth at until the modes enter the Hubble radius. Thereafter the modes decay as normal vectors, 
with some gentle oscillatory features. 



spectrum tends to affect the amplitude of second-order modes on large scales at the level of a few percent [IJ, |20j ■ 
Viewing x as time for constant fc, we see that the modes grow as t], peak when inside the Hubble radius and decay as 
?7~^. While the modes are decaying there are faint oscillations as shown in the top panel of the Figure. This Figure 
can also be interpreted as the power spectrum at fixed time, showing the usual features. 

It is worth mentioning that we have taken the upper limit of the k' integral to be infinity. In reality there is a 
cutoff from the end of inflation, at = 77* , corresponding to modes which are inside the Hubble radius at that time, 
/c* = l/?7*, so giving a finite upper limit to the w-integral, = k^/k (although this is on very small scales in reality). 
This causes a break from linear scaling in x in the power spectrum for x < and we may analytically find the 

leading behavior of Eq. ((^ is Vv ^ ^v^x^ . Why is this the case? 



In the interacting delta function case we saw that modes are efficiently produced - and grow like - when both 
modes are outside the Hubble radius; once one is inside and the other outside there is effectively no generation 
of vectors. In power-law case, then, the modes which are generating vectors are those outside the Hubble radius, 
providing an effective cutoff to the w-integral of f ~ so giving us growth oc when 3> l/x. When we have 
the cutoff v.if on the other hand, for early times when 77 < l/fc.^, all relevant modes are outside the Hubble radius, 
and interact coherently giving us growth oc 77^. For r\ > 1/fc*, modes which have entered the Hubble radius no longer 
contribute to the generation of modes outside the Hubble radius giving weaker growth cx x. Of course, we are not in a 
position here to analyse times before inflation ends, but this helps us understand why we have the x-scaling behaviour 
we do. 



An important feature of inflation is the lack of vector modes: if they were observed to have a similar spectrum 
and amplitude to the scalars then this could prove difficult for inflation, and lend favor to other theories of the early 
universe, e.g. Pre Big Bang scenarios and Ekpyrotic models [21I, [5l|, [s^l ■ However, there is a x^-distribution of vectors 
produced by inflation as a consequence of the non-linear interaction of scalar modes, which has received relatively 
little attention to date. 

We have investigated the generation of vector modes induced by primordial density perturbations during the 
radiation dominated era. Performing a perturbative expansion to second order, we isolated the scalar terms which 
source the vector perturbations. We then calculated the power spectrum of the metric vector mode, and analysed 
its form. In order to understand the generation of modes we investigated individual scalar modes generating vectors, 
and demonstrated that, contrary to the case of gravitational waves, vector modes can't be generated by an isotropic 
distribution of scalars of a single wavelength, owing to the spin-1 nature of vector modes: rotational degrees of freedom 
must be generated by scattering of non-parallel input modes. We then demonstrated that vectors are generated by 
modes of differing wavelength whenever one of the two scalar modes is entering the Hubble radius. The amplitude 
of the generated modes depends on the ratio of input wavenumbers; maximum power is generated when the modes 
are not too widely separated. After investigating the generation of modes, we then presented the power spectrum for 
scale- invariant scalar modes, displaying our results in terms of the variable x = kr]: i.e., they may be interpreted the 
temporal evolution of a single scalar mode, or the power at a flxed time. Interestingly the maximum power generated 
is the same at all times, but the position of this peak changes with wavelength, such that x ^ 1 ^ k ^ I/77. This 
is due to the fact that the modes are efficiently generated as they enter the Hubble radius, and aren't generated 
significantly while outside. 

There are some open questions raised by the study presented here. In particular, it is not clear how the power 
spectrum for Si we have calculated will be related to observable quantities. It is gauge invariant so must be observable, 
by the results presented in [sl]; it also represents all possible degrees of freedom of vectors generated by scalars, under 
the conditions laid out here. Thus, although there may be a 'better' variable, it must be related to Si by quadrature 
(plus some further scalar-squared contributions). The effects we have presented here will have interesting implications 
for a variety of phenomena such as the CMB; the issue of how significant is left for future work. 
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In this section we consider the extraction operator discussed in this paper. We start by defining the Fourier basis 
used for the purposes of harmonic decompositions. An arbitrary scalar in real space can be expressed as a Fourier 
integral 
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APPENDIX A: THE EXTRACTION OPERATOR 
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A divergence free vector in real space can then be expressed as a Fourier integral 
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where ea{k) and ea{k) are orthogonal parity vectors, which are also orthogonal to k. Similarly, a transverse traceless 
tensor (a tensor mode) in real space can be expressed as 



Tab{x,v) = I d'k [T{k,7j)qab{k)+f{k,lj)qab{k)] 6*^ 
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where the two polarization tensors qat and qab can also be expressed in terms of the parity vectors Ca and Ca, and are 
orthogonal to A;. An arbitrary symmetric trace- free spatial tensor in real space 
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which has explicit scalar, vector and tensor contributions. We can also express Aat as a Fourier integral 
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The even-parity vector contribution can then be extracted in Fourier space by applying the following operator 
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Applying the operator Vi"-'' to Eq. IA4I gives 

A\''\x,7j) = Vi-'Aab 

2i 
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d^x' k'- [ei{k')e\k') + ei{k')e'{k')] e''''<^-^'^ Aab{x' (A6) 



The extraction of the vector component is made possible by taking advantage of the various properties of both the 
parity vectors (ca and Ea) and the polarization tensors {qab and qab)- We now consider as an example one possible 
contribution from the first order squared terms. For simplicity we start with contributions from a term made up of the 
product of first order scalars and do not reconstruct in real space, looking only at the Fourier amplitudes. Consider 
a term of the type 



1 



(27r)3 



The Fourier amplitude of the vector part of this is then 

2i 



J d^ki 77)6*'=-=^ j d^k2 ^k2,v) {~k2ak2b) e"^" 

k- [e'' (fc) + (fe)] k2ak2b 



(27r)3/2 

where we have carried out a real space integral and a /c-space integral. 



fc2 



(A7) 



(A8) 
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